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PREFACE ' 

♦ ' 

The report consists of two parts i Part X Is a study of lumped parameter 
models and their ability to predict the stiffness and damping characteristics 
of gas film bearings and dampers. This effort* was assisted by Shun-Lung Chao 
and represents his Master of Science thesis. Part II is an experimental 
investigation of the separation bubble near the inlet feeding region of 
inherently compensated gas bearings. Contributions to this work include 
Don Dye, a graduate student and Keith Logue, an undergraduate student in 
Mechanical Engineering. 
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PART I: 

A LUMPED PARAMETER MODEL FOR 
GAS FILM BEARINGS 



iv 


.ABSTRACT 

A lumped parameter model is developed to determine the stiffness and 
damping characteristics of inherently compensated gas film hearings. The 
model relies on the average static pressure over a one-dimensional strip 
bearing. Results of the model are compared with known computer solutions 
for the distributed strip and a two-dimensional square bearing. The 
results for the stiffness agree well with the computer solutions although 
the model proved to be inadequate for predicting the film damping. 
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NOMENCLATURE 


SYMBOL 

DESCRIPTION 

* 

Denotes actual variable 

0 

Denotes average value 

Al 

Area of the central region 

\ 

Area of the bearing pad 

C* 

Damping 

C 

Dimensionless damping 


Orifice discharge coefficient 

D 

d/dt* 

So 

Gravitational acceleration 

h* 

Film thickness 

h * 

Average film thickness 

h 

it 

Dimensionless film thickness ( h / 

k 

Ratio of Specific heats (Cp/C^) 


Stiffness 

Ks 

Dimensionless Stiffness 

L* 

Bearing width 


Mass flow "to -the -central .region , • 

Mil* 

Mass flow through sill region 

M/ 

Hass flow through the inlets ■ • 

Min* 

Mass flow into the bearing pad 

Mout* 

Mass flow out of the bearing pad 

“’o* 

Average mass flow 



X 


SYMBOL ■ 


P* 



^ro 


r 


R 

Sd* 


t 


*. 


t 


T 



W* 


W 

W 



NOJtENCLATURE (continued) 

DESCRIPTION 

Pressure 

Ambient pressure 

Dimensionless pressure ( P' /P^ ) 

Dimensionless pressure downstream of inlet 
Dimensionless average pressure doTOstream of inlet 
Supply pressure 

Average pressure across the sill region 
Inlet span to total bearing span ratio 
Gas constant 

Dynamic stiffness ( ~ AW / Ah”) 

Time 

Dimensionless time ( t ) 

Gas temperature 

Film volume of the central region 
Film volume of the sill region 
Bearing width 

Total bearing load capacity 
Dimensionless bearing load 
Average bearing load 
Cartesian coordinate 
Gas density 

Film density of the. central region : 

t • 



SYMBOL 

/?i* 

(T 

£ 1 * 

A 
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'NOIILNCLATURE (continued) 

DESCRIPTION 

Film density of the sill region 
Gas viscosity 
Squeeze Number 
Excitation frequency 
Restrictor coefficient 

Ratio of the central area to the bearing area (Aj/Ap) 
Amplitude of disturbance 



• CHAPTER I 


, INTRODUCTION • . , 

t • ’ , ' ' ’ ‘ 

The range of application for externally pressurized gas lubricated 
bearings has been increasing steadily in recent years. The reason for 
this trend is the many advantages of gas lubrication. Gas bearings 
offer extremely low friction due to the low viscosity of the gas; so 
mechanical efficiency is high, but heat generation is low, especially 
in low or medium speed rotation. Because of the absence of metal-to- 
inetal contact which promises freedom from wear, the life of the bearing 
is long. Also, gas bearings provide silent and vibrationless running 
and are able to operate in extremely high or low temperature. In 
addition, because of the distributing action of the gas film, gas 
bearings can reduce the effect of local errors of geometry and produce 
an accurate mode of motion. According to these advantages, gas 
bearings have been successfully applied in increasing amounts in many 
industrial applications over the past few years. 

Although gas lubrication offered many advantages, two major 
problems exist:. 

i 

(1) , the existence of a self-excited film whirl instability [l}* 

in which machines can not be run in excess of the whirl onset 
speed; 

(2) the limitation of the bearing load capacity which usually 
results in severe damage to the bearing surfaces by overload- 
ing, . 


.^Numbers in brackets refer to REFERENCES. 
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The response of gas bearings la affected by the dynamic 
characteristics of the lubricating film. To study this problem, the 
gas film can be approximated as a linear spring-dashpot system . 
subjected to a small amplitude, periodic load disturbance. Thus, 
stiffness and damping are important parameters describing the djTiamic 
response of this system. In order to determine the- dynamic bearing 
load, the Reynolds' equation must be, solved for the pressure distribu- 
tion" in the gas film. The . Reynolds ' equation is a nonlinear, time- 
dependent, parabolic, partial differential equation; so it is difficult 
to find exact solutions for this equation.. For obtaining approximate 
solutions, small-perturbation theory is often combined with finite- 

f' ' . ■ U . 

difference methods; Limited .use has been made of so "called " lumped 
parameter" methods, but their accuracy is not clearly understood. A 
review of the literature discussing the above methods will now be given. 

1.1 Review of literature 

Richardson [ 2 ) analyzed the static and dynamic performance of gas 
bearings based on lumped parameter methods. He derived the relation- 
ships of the bearing parameters to be used for design Information. 

Licht and Elrod [ 3 ) studied the- stability of externally pressurized gas 
bearings, using distributed .parameter methods. Their results when 
compared with those obtained by Richardson 'using lumped parameter 
methods showed a marked divergence among the limiting values of 
parameters which affect the stability of the bearing. Stiffler[4j 
used distributed parameter methods to study an inherently compensated, 
multiple-inlet, circular thrust bearing. In his analysis, the 
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Reynolds' equation was solved by small-perturbation theory, and the 
results were given for the bearing coefficients (stiffness and damping) 
as functions of inlet location, excitation frequency, supply pressure, 
and restrictor coefficient. Mullan and RichardsonfS) used ‘lumped and 
distributed parameter methods to develop the dynamic characteristics 
of ari inherently compensated gas journal bearing and compared the 
analysis with experiment. 

Lund [6) carried out a theoretical analysis of the load-carrying - 
capacity of the hydrostatic gas journal bearing vrith jouraal rotation 
and vibration,' Results were obtained from a first-order perturbation 
solution by assuming small eccentricity ratio and small vibration 
amplitude, Also,^ Lund (?) analyzed the tKreshhold of instability for 
a rigid rotor supported in pressurized gas journal bearing based on . 
a first-order pertubation with- respect to the eccentricity ratio and 
makes use of the "linearized Ph" methods. Numerical results are given 
as a function of supply pressure ratio, feeding parameter, and eccen- 
tricity ratio. > 

Stiffler and Smith [s) analyzed the dynamic characteristics of- 
an inherently compensated,, square gas film bearing using small pertur- 
bations of the Reynolds* equation. They obtained the results by 
finite-difference methods and presented the design curves for the load 
capacity, mass flow, stiffness and damping as a function of squeeze 
number, supply pressure, restrictor coefficient, and inlet location. 
Stiffler and Tapia (9J presented a theoretical analysis of the 
amplitude effects on the dynamic characteristics of an inherently, 
compensated, infinltively long strip, gas thrust bearing. 



A 


Squeeze number (dimensionless frequency) is an important parameter 
in understanding the dynamic response. Sadd and Stiffler (lO) 
analyzed gaseous squeeze film dampers at low squeeze numbers and 
determined the effect of periodic disturbance amplitude on the 
dynamic performance. The results are useful in predicting the 
behavior’ of externally pressurized gas bearings uiider certain 
conditions. Earlier, working with circular squeeze films, Salbu Cn3 
presented computer solutions of the pressure distribution as a 
function of the squeeze number. A superambient average film pressure 
was introduced by compressibility effects at high squeeze numbers, so 
that a gas film can support a load while operating on squeeze effects 
only. These results were supported by experiment. 

■ ' t ' , 

' * ‘ " 

1.2 Statement of the Problem 

The general problem of obtaining the dynamic characteristics of 
externally pressurized, inherently compensated, gas thrust bearings 
is the solution of the Reynolds’ equation, a non-linear partial 
differential equation in space and time. Solutions of this equation 
by finite-difference techniques have presented difficulties in the 
past, particularly related to convergence time. Furthermore, the 
multitude of parameters have limited the amount of information that 
can be received from the computer and illustrated for the designer. 
The purpose of this thesis is to develop a lumped parameter model for 
these bearings in analytical form. Such a model can provide a quick, 
though approximate solution for the designer as well as provide a 
check on the more accurate finite-difference solution if it is so 
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desired . 

The following chapter on analytical methods contains txro sections 
: (1) a finite-difference solution of the strip bearing, Figure 1, 
which has been condensed from REFERENCE [9J-; (2) a general lumped 
parameter model for all bearings. The solution of the strip is 
outlined here for several reasons: (1) parallels can be observed 
between the ti70 methods; (2) many of the basic equations are pertinent 
to both methods, including dimensionless parameters; (3) Itimped 
solutions of more complicated tiro-dimensional bearings can be inferred 
from a model of the strip which includes aspects of the exact solution; 
(4) comparisons between two solutions can give some insight into the 
accuracy of the lumped model. 








CHAPTER II 


ANALYSIS OF AN INHERENTLY COMPENSATED STRIP THRUST BEARING 

2»1 The Strip Bearing (Condensed from RlilFERENCE (9)) 

To develop the theoretical analysis of the strip bearing, it is 
divided into two different regions (Figure 1) ; the central region or 
region I, and the sill region or region >11. In order to simplify* 

, i ' f , ^ ^ ^ 

the problem, only one-half of the strip bearing is considered 
because the geometry is symmetric about' the centerline. Further— 
more, the inlets are assumed to be sufficient in number so that • 
they may be considered as an inherently equivalent line source. It 
is also assvmied that the Length of the strip bearing in the Z-axis 
is infinite. 


2,1.1. Reynolds* Equation 

Before solving for the bearing mass "flow, bearing load 
capacity, bearing stiffness, and bearing damping, the pressure 
distribution across the bearing first has to be determined. The 
pressure distribution in the clearance of an externally pressurized 
bearing is governed by the Reynolds’ equation [iz) ; 


( h*' II ) + ^ ( h*) (2- 

Of:' 


* 


■9t' 


la) 


^P 

6x* 3x* 9z* 

where the starred symbols denote real varxables. Since two sur- 
faces of the bearing are parallel and rigid, the film thickness, - 


h*, is uniform across the bearing. 

For the strip bearing, the pressure distribution is- described by 
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the reduced Reynolds* equation: 

A ' "* V’ ^ ( /=* h*) ' , (2-ib) 

ox o X <Jt 

In order to simplify the analysis, the film is assumed to behave as 
an ideal gas with constant specific heats. It is also assumed that 
the flow is isothermal x^ith P*//3*=constant. 

To solve equation (2-la), the following nondimens ional variables 
are introduced f 

X = X* / . (2-2a) 

h =* h*/hQ (2-2b) 

P = P*/Pa*' ' , ■ (2-2c) 

, t =» t*ri*. 

■ ’ ' ‘ ‘ 

Thus, equation (2- lb) can be written as 


>2 


>2^ _ 2(T a 


( P ) - 
3x2 ' h^ 9t 


(Ph) 


(2-3) 


where the squeeze number, (p, is defined as 


(T = 




(2-4) 


and the subscript symbol,©, denotes average condition. The film 
thickness, h, is a function of time because it depends on the load 
disturbance. 

To solve equation (2-3) , the boundary conditions must be knoxra. 


The average pressure doxmstream of inlets is given by Pro » ^nd the 
pressure at the edges of the sill region is ambient. Thus, the 
boundary conditions, x^hich are normalized by the ambient pressure 



are 

P = Pj,Q at X‘s» a (2-5a) 

P « 1 at X = 1 - (2-5b) 

' , ' ^ I 

Applying, equations (2-5a) and (2-5b) to equation (2~3), the 
average pressure distribution .through the bearing can be determined as 


P « P 


ro 


0< a 


(2-6a) 


p= 1+^.^jcali. j (1-x) 




x<l 


(2-6b) 


where 'a' is the dimensionless distance from the centerline to the 

• ' , i 

inlets. ‘ ■ 

^ ‘ ‘ ‘ ’ ' ' M ‘ ' 

2. ,1.2, Mass Flow > v • ’ ’ 

From the continuity condition at the inlets, the following rela- 
tion has to be met: 


Mj.* + « Mu* 


(2-7) 


where 


Mj.* = mass floX7 through the inlets 

Mj* ~ mass flox7 through the central region 

Mjj*=* mass flow. through the sill region 

The mass flov; across two parallel, long strips is given by 


Cons tan tlnescu (12) as 
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fL* * ap* ^ * 

rP 12,11 Sx*' 


(2-8) 


Since, the flow; is isothermal, equation (2-8) becomes 


. = - 


3 ■ 2 
L*h* Pa* 


. 2^RTw* . dx 


(pZ) 


(2-9) 


Thus, the average ‘mass flow through the central region and the ■ 
sill region can be determined from equation (2-9) by applying equations 
(2-6a) and (2-6b) to obtain 


■ Mj* = 0 


and 


« L*h*^Pa*^ (Pr^ ~ D 

1^RTt7*- (1 - a) 


( 2 - 10 ) 


( 2 - 11 ) 


The mass flow through the inlets (l3l is described by 


M * 
r 


2CnL*h*PaPs 

(RT)i2 ^ 






Ps + 1^ 


(2-12) 


or 
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Since tlie average mass flox^ through the inlet and the sill are 
the same, the average pressure downstream of the orifice, Pj.q, Is 
established by 



ro 



( 2 - 14 ) 


or 




k 


( 2 

ik+li 


( 2 - 15 ) 
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where the restictor coefficient,/^ 


is defined by 


■ ' ' ^ ' ' 24CpyHW*(l-a) /2gpkRT\ ^ 

i ' VaX*^ : k-1 / ( • ‘ (2-16) 

■ ■ ’ ■ ' ' • , I \ 

' » i ^ 

The restrictor coefflcient,/\ils a nondimensional parameter which 
indicates the ratio of the resistance of the mass flow through the 
sill region to the mass flow resistance across the inlets. For 
critical flow, Pj-^ is obtained directly from equation (2-15), But 
for subcriti'cal flow, P^-q is solved by the. Newton-Raphson method. (see 
APPENDIX for tables relating Pj-q to A) since equation (2-14) is 
implicit for Pjo. After Pj-^ is knox-m, the pressure distribution, 
equations (2-6a) and (2-6b), can be determined. 


2,1.3 Bearing Load Capacity 

The bearing load capacity is given by 


W* 


fw* 


(P* - Pa*) L*dx* 


2w*L*P 


'•p^*r (P-l)dX 


(2-17) 


where the pressure P is- attained' from the reduced' Reynolds ' equation 
(2-3), Thus, the dimensionless load capacity is derived from 
equation (2-17) as 


W* 


W = 


2w*L*?a* 


o 


( P - 1 ) dx 


(2-18) 
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Substituting equations (2-6a) and (2~6b) into equation (2-18), ■ . 
the average load capacity is obtained as 

r / 3 \ ' 

¥ « a(Pr - 1) + i 1 (1 “ a) (2-19) 

3 \P^ - 1/ / 


Thus, the average pressure across the sill region can be found 
from equation (2-19) as 



(2-20) 


2.2 Lumped Parameter IJbdel 

Since the gas film behaves as a simple, linear, first order 
■system for small disturbance amplitudes, the load changes, vrhich are 
usually periodic in nature, exhibit components in phase and ninety 

degrees out of phase with displacement. Thus, a dynamic stiffness is 
defined as 

Sd* = - 

Ah* 

(2-21) 

The bearing load capacity is a function of the pressure down- 
stream of the inlets only, so 


AW* 

AP * Ah* 
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•AU* 

APj-* 


=• 2w*L* 




dP, 


( 2 - 22 ) 


where the bearing load capacity is given by equation (2-19). 

■ The conservation of mass flow requires that the net mass flow 
into the bearing pad is equal to the change of mass stored in the 
bearing pad. Thus, 


M. * ■- tl * 
xn - ' out 


(mass) 


dt* 


pad 


d d 

— — (/0*V ) + 

dt* ' central region dt* 


(jO*V) 


sill region 




(2-23) 


where 


D 


d 

dt* 



Vli 


■'= film density of the central region 
= film volume of the central region 
*= film density of the sill region 
= film volume of the sill region 
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Linearizing equation (2.-23) for small changes from the steady- 
state conditions, equation (2-23) becomes 

“(Vi) + '>10 + ''iio“‘/’n^l 

(2-24) 

where the symbol, o, -denotes average condition. Also, the following 
relations are introduced: • 


P * -J-4 P* 
^ro 


ho* +Ah* 


A 


Ai/Ap 


(2-25a) 

(2-25b) 

(2-25c) 


where Aj and Ap are the area of the central region and the total 
bearing pad respectively. 

Since the average .pressure across the sill region is found from 
equation (2-20) , the following relations are derived. 




lo 


p * 

= (-£a^) p 
RT 


ro 


D(Vi) = A Ap D(4 h*) 


^lo ” ^ ^p^o 


D 


fl*') - (- 5 - ) D (4 P/) 




llo 


) (' 


p 3-1 

^ro 


^o - 1 


(2-26a) 
^ (2-26b) 

(2-26c) 
(2-26d) 


RT 


(2-26e) 



16 


B(vn) 

• (1-X) Ap B(4h*) 

(2-26f) 

'^IIo 

- a-X) Aphp* 

(2-26g) 



2 1 

“( ) 

3 RT 




(Pro^, - 1)' 


. ^ . 

D (AP„ ) 


(2-26h) 


The mass flow into and out of the bearing pad is a function of 
film thickness and the pressure downstream' of the inlets. Thus, the 
following relations can be found: 




9^in \ * 

|Ah + 


a h' 




A / r 


(2-27a) 


4M 


^9 'A 




out 


but • 


.3 h' 




(2-2 7b) 


Substituting equation (2-26) and (2-27) into equation (2-24) , the- 
following relation is obtained: 


A Pi 


Ah 


V. 


A r 


RT 


3 '• 


1 ? “'-1 


P ^-1 
^ro ^ 


D + 




n 


v-ah 


9 h 


A 


A h. 


RT 


3 


' Pro®ro^-3Pro+2) 


(Pro^-l)' 


3Ml * 


D + 


‘xn 


s^app* apj 


(2-28) 
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Equations (2-22) and (2*.28) are introduced into equation (2-21), 

J 

then 

• * 

®d 

where 





; 

A linear approximation of the bearing motion is 
AW* » ” C*D[A-h*) 


(2-31) ■ 


(2-32) 


(2-33) 
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- Thus, 



■ A - 

= •+ C ,.D 

3 • i 


(2^34) 


V A 

where K and C represent the stiffness -and- damping of the bearing 

s ' r ^ * 

respectively. 

For a periodic load disturbance, the D operator in equations 

• A 

(2-29) and (2-34) can be replaced by j^operator. Then, 


Sd = k 


a+n* riT2) + (?i “'^2) 


* 




(2-35) 


Sd*' = k/+ c*jif 


(2-36) 


From equations (2-35) and (2-36) , the stiffness and damping can 
be determined as 


, ks* ( j +n *^T{r2) 

■ ^ (1 (2-37) 

i T2^ . (2-38) 

For certain conditions, equations (2-37) and (2-38) reduce to 
the following:* 

■ t.tl'an ^ ^ ^ (2-39) 
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K- 


Tl k, 

T 2 


when 


Q. » 


X 


and 




(2-40) 


kg* (Ti -r2> 


A 1 

when ^ « -=5“ 

Co 


(2-41) 


The nondimensional stiffness and damping are defined by [9^ 


K 


K * h* 
_£L Q- 


2w*L*P^*(Pg - 1) 


(2-42a) 


k * h * 

Itg IIq 


2w*L*Pg*(P^ - 1) 


(2-42b.) 


(t^)' 

-V 


(2-43a) 


(2-43b) 


From equations (2-10) , (2-11), (2-13), and (2-19), the following 
relations can be obtained: . 


dW 

dP 




= 2w*L* 


ro 


2(l-a) /P 


a + 


ro 


3 


3P 

~’^ro 


2(P ^ - 1) 

ro 

(P ^ - 1) 

' ro ' 




(2-44) 
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out 


8 h 


K * p 

“o » ro 



^^Wt' 




= 2 


h P 

“n 9 ^ 


ro 


^ro’^o 


Pa*(?ro^ 


- 1 ) 


where ^ L*Pa*^ho* 




✓«RT w* 


™o 


(2-45) • 


(2-46) 


(2-47) 


(P ^ - 1) 
ro 

m = 

•12(1 - a) 


and 



9 h" 


^ro 



(2-48) 


(2-49) 







or 


a%n 

3h* 


in. 


ho*,P 


ro 


(2-51) 
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Thus, after the inlet location, supply pressure, average pressure 
downstream of the inlets are knoim, the dimensionless stiffness and 
damping can be found from equations (2-53), (2-54), (2-55), and (2-56). 



CHAPTER III 


RESULTS AND MODEL COMPARISONS POR THE STRIP BEARING 

^ \ 

The stiffness and damping are obtained in this study as functions 
of restrictor coefficient and the supply pressure with inlet location 
(a-0.5' in this study) fixed, .They are compared with computer solutions 
by Stiffler and Tapia [9] at low amplitude of disturbance (£=0.1) and 
small squeeze number ((T=0,1). The amplitude effect on the stiffness 

. * I < , * - 

and, damping can be neglected at low amplitude of disturbance. The ; 
stiffness and damping are insensitive to small squeeze numbers ((T<'.4) , 
but the stiffness increases and the damping decreases x^hen the squeeze 
number exceeds this value [9 J . Design curves for the stiffness and 
damping are discussed and compared beiox'? for the lumped .parameter 
approach and the computer solutions. 

3.1 -Stiffness 

The dimensionless stiffness is a function of the restrictor 
coefficient, A , and supply pressure as shorn in the Figures 2-6. The 
stiffness is very sensitive to the restrictor coefficient. When the 
restrictor coefficient becomes small, the stiffness approaches zero. 

The stiffness ‘ reaches the maximum in the range 1.3</\<2.2 where the - 
mass flow through the inlets is critical. As the restrictor coeffi- 
cient becomes large, the stiffness approaches zero again. 

If the stiffness is compared with computer solutions by Stiffler 
aixd Tapia (9J as shown in Figures2-6, it can be observed that the 
results of the lumped model agree well x^ith computer solutions. 



RESTRICTOR COEFFICIENT 


Figure 2. , Strip Bearing-Dimensionless Stiffness versus Restrictor Coeff icien 
(Pg=1.5-,(r=0.1, a=0.5) 
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Figure 3, Strip BeariTig-Dimensionless Stiffness versus Restrictor Coefficient 
(Pg=2,(r=?0.1, a=0.5) 
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Figure 5* Strip Bearing-Dimensionless Stiffness 
(P3-6,(T=0.1, a-0.5) 
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Figure 6. Strip Bearing-Dimensionless Stiffness versus Restrictor coefficient 
(Pg=10,(T=0.1, a=0,5) • 
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3.2 ■ Damping . 

The relationship between the diiaens ionless damping, C, and the 
restrictor coefficient, A » is shown in Figures 7-11. It can be obser- 
ved that the dimensionless damping is highly dependent on the 
restrictor coefficient and the supply pressure. ■ As the restrictor 
coefficient approaches zero, the damping approaches a constant. As the 
restrictor coefficient approaches infinity, the damping approaches a 
smaller constant. These values represent the extremes for pure squeeze 
films ,[10)' i.e. no inlets (A=0) or Inlets open to the ambient (A=°o). 

Comparing the damping between computer solutions by Stiffler and 
Tapia(9} and lumped parameter approach by this study, it is observed 
that the lumped model agrees well ^d.th computer solutions for the range 
of lower restrictor coefficients,/i<1.0, but differences occur in the 
range of larger restrictor coefficients. These results indicate that 
a lumped model for the damping is not' entirely adequate. Sadd and 
Stiffler (10} have shown that pure squeeze films (A =0,oo) display quite 
different damping depending on the geometry. However, the lumped 
model would display no difference between a strip and a square, for 
example, if their corresponding film volumes were identical. 

3.3 Squeeze Number 

To introduce the .squeeze number, c, into equations (2-37), 

(2-38) , (2-42a) , and (2-43a) , the dimensionless stiffness and damping 

can be expressed as follov7s: 
when 

■ 

Pg Vk + I/ 
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Figure 8. Strip Bearing-Dimensionless Stiffness versus Restrictor Coefficient 
- (P3«2,<r=0a, a =0.5.) 
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Figure 9. Strip Bearing-Dimensionless Stiffness versus Restrictor Coefficient 
(Ps=^»<S'“0.1, 1=0; 5) 
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Figure 11 • Strip Bearing-Dimensionless Stiffness versus Restrictor Coefficient 
(Ps=10,cr=0.1, a-0.5) 
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where 



^ The effect of squeeze numbers on the stiffness and damping is 
shown in Figures 12-17. It is found that low squeeze number (<T<1) has 
little effect on the stiffness and damping. Hpwever, the stiffness 
increases and the damping decreases for large squeeze numbers. 

Computer solutions are available for Ps=10,A“l and are presented 
in Figures 12,13 also. The break frequency for the lumped parameter 
damping in Figure 15 appears to be realistic in spite of the short- • 
comings reported in Section 3.2. 



CHAPTER IV 


ANALYSIS OF THE SQUARE BEARING 


4,1 The Square Bearing (Condensed from REFERENCE [s]) 

To analyze the square bearing, it is also divided' into two 
regions (Figure 18) as the previous strip bearing. From the symmetrj' 
of the geometry, one-quarter -of the square bearing is considered here, 
and an inherently equivalent line source is used to replace the 


discrete inlets. 


V 


( ► 




The Reynolds’*' ’equation, equation (2-1), can be normalized using 
previous dimensionless variables by using the length L* in place of 
width w . • ■ 

Thus, ' • 


B <■>') + B (p^) - ^ - 

3z^ h3 ut 


(Ph) 


(4-1) 


where the squeeze number, (T, is defined as 
* 

L2/iSl L* 


(T = 


h_*^ * 

I’o 


W-2) 


The average mass flow across the sill area is given by REFERENCE 


( 8 ) 


as 


* . *8 


m. 


* Pg 




2«RT 


where 


(Pro^ - 1) F 




(4-3) 


(4-4) 


and F is given for the three different span ratio in the following 


table: 
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F 


0.4 

0 . 6 . 

0 . 8 . 


1.83 

3.48 

8.44 


^ y From the flow continuity, equation (^2-7) , for the average mass 

• ' i ' ; . ; ’ 

flow across the bearing, *the> ‘average pressure do\«istream of the 

I 

orifice, Pro» is determined by equations (2-14) and (2r-15) except that 
now the coefficient, A » Is 


6Cp NTi: d/^ A2gpkRT\ ^ . 

^ ■ ?s Pa* F U - 1 ) 


(4-5) 


The bearing load capacity is given by 

fL*/2 fL*/2 . * ^ 4. 

W* - 4\ \ ( P - P. ) dx* dz* 

d 


« -4L*^ Pa* (P - 1) 

Thus, the dimensionless load capacity is 


W = 

L* P * 
a 


(hfh 

= 4) \ ( P - 

Jo Jo 


1) dx d z- 


(4-6) 


(4-7) 


In order to apply the lumped parameter model to the square 
bearing, it is assumed that the average pressure on the', sill region 
is closely approximated by the average sill pressure exhibited by the 
strip. Then the average load- capacity is expressed by - 
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where the area ratio, Is given in terms of the span ratio. Figure 
18, as ■ 


A “ 


(4-9) 


The previously developed equations for the strip can be used for 
the square bearing providing one interprets the results in terms of 
equation (4-4), (4-5) and (4-9) « 



CHAPTER V 


• RESULTS AND MODEL COMPARISONS FOR TH'e SQUARE BEARING ‘ 

The stiffness and damping are obtained in this study as functions 
.of restrictor coefficient and supply pressure vith inlet location 
( r = 0,6) fixed. They are compared with computer solutions by 
Stiffler and Sm.ith [8) at small squeeze number (<r® 0’. 1) since the 
stiffness and damping are insensitive to small squeeze numbers. Design 
curves for the stiffness and darapling are discussed and compared below 
for the lumped parameter approach and the computer solutions. . 

5.1 Stiffness 

If the stiffness is compared with computer solutions by Stiffler 
and Smith [Sj as shorn in Figures 19,20, it can be observed that the 
results of the lumped model agree well with computer solutions. 

However, the sti ffnes s for Pg=2 is greater than it for Pg=6 in’ the 
range 1.5<A< 5.5. This fact indicates that the lumped parameter model 
can not exactly describe the performance of the square bearing, but 

the agreement is snfficient.fo;: most engineering purposes. 

* * } 

f 5*2 Damping 

^ V ' 

Comparing the damping between computer solutions by Stiffler and 
Smith [8] and the lumped parameter approach by this study, the results 
show that the lumped parameter model is not adequate enough to 
describe the performance of the square bearing. Sadd and Stiffler 
(l0]have show that the damping for the squeeze bearing is approxi- 
mately one-half of the damping for the strip bearing in the case of 
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Figure 22. Square Bearing-Dimensionless Damping versus Restrictor Coefficient 
(Ps=6,<r=0.1, r=0.6) 
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pure squeeze film (A ~0) . So if the damping from the lumped parameter 
approach, is adjusted by the factor 1/2, the results are improved but 
the agreement could not be described as sufficiently accurate. 

5.3 Squeeze Number 

Comparing with computer solutions[S } in Figures 23-26, it is 
observed that the break frequencies for the lumped parameter model 
appear to agree well x^ith the computer model in spite of the inadequacy 
of the model as indicated in Section 3.2. According to this result and 
observing equation (2—38), the term which indicates the compressi- 
bility effect of the bearing film alone controls the break frequency. 
Thus, the terra which indicates the volume effect of bearing film 
xTOuld appear to be the target variable in attemps to improve the lumped 
model for damping. 
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CIIAPTER VI- 


CONCLUSION 

The dimensionless design parameters. stiffness and damping obtained 
from the study of the strip bearing by lumped parameter methods are 
compared with those of Stiffler and Tapia [9 ) who used finite- 
difference methods. It is observed that the stiffness obtained from, 
this study agrees well with computer solutions; however, the damping 
works well for the range of low restrictor coefficients only, A510, 

In discussing the effect of squeeze numbers on the stiffness and 
damping, the break frequencies for the lumped parameter damping appear 
to be realistic in spite of the inadequacy of the model. 

The stiffness and damping obtained from the study of the square 
bearing by lumped parameter methods are compared with those of 
Stiffler and Smith [8] who used finite-difference methods. The 
stiffness obtained from this study is adequate when compared with 
computer solutions,; however, the limped model damping results offer 
poor agreement. If the damping from lumped. param.eter solution is 
adjusted according to Sadd and Stiffler’s work on pure squeeze films 
(lO), the results are improved but remain Insufficient for predicting 
performance. On the study of the effect of squeeze number on the 
stiffness and damping, the break frequencies can be predicted well in 
spite of the inadequacy of the lumped model. 

Since the lumped parameter model is an approximate technique to 
study the performance of the bearing, it is not surprising that the 
results differ from the computer solution. The corresponding film 
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volumes of the strip and square bearing are comparable in this study; 
so the ♦ lumped model displays no difference betueen a strip and a 
square^ contrary to known results for pure squeeze films. This is 
one of the reason why much difference exists for the damping of the 
square bearing compared with computer solution. Because the predic- 
tion of the break frequencies for the lumped damping agrees well with 
computer solution, it can be concluded that the compressibility effect 
of the bearing film is well described by the lumped model, but the 
volume effect of the bearing film directly affects the accuracy of the 
lumped model. This suggests a starting point in future attempts to 
improve the lumped models for predicting bearing performance* 

Regardless, the lumped model predicts the stiffness performances 
of -both -the strip and- square bearing. Stiffness remains the single 
most important parameter in the design of rotor-bearing systems. 



APPEISDIX 


TABLES FOR DIMENSIONLESS AVERAGE PRESSURE 
DOWNSTREAM OF INLET Pj.q VERSUS RESTRICTOR 


COEFFICIENT A 
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TABLE 1 


• P - versus A for P » 1.5 
ro s 


A 

Pro 

A 

Pro ■ 

A 

Pro 

30.0 

1.49821271 

3.0 

1.39086568 

1.1 • 

1.22797903 

25.0 

1.49743451 

2.9 

1.38608333 

1.0 

1.21224958 

20.0 

1.49601475 

2.8 

1.38102095 

0.9 

1.19626465 

15.0 

1.49300268 

2.7 

1.37565945 

0.8 

1.17926772 

13.0 

■ 1.49076996 

2.6 

1.36997833 

0.7 

1.16120100 

10.0 

1.48478685 

2.5 

1.36395567 

0.6 

1.14200700 

9.5 

. 1.48325201 

2.4 

1.35756799 

0.5 

1.13628988 

9.0 

1.48147882 

2.3 

1 . 35079020 . 

0.4 

1.11037100 

8.5 

1.47941673 

2.2 

1.34359550 

0.3 . 

1.08383246 

8.0 

1 . 47700148 ' 

2.1 

1.33595521 

0.2 

1.05662759 

7.5 

1.47415061 

2.0 

1.32783875 

0.1 

1.02870352 

7.0 

1.47075716 

. 1.9 

1.31921364 

0.09 

1.02586931 

6.5 

1.46668082 

1.8 

1.31004522 

0.08 

1.02302725 

6.0 

1.46173548 

1.7 

1.30029668 

0.07 

1.02017727 

5.5 

1.45567120 

1.6 

1.28992908 

0.06 • 

1.01731930 

5.0 

1.44814827 

1.5 

1.27890107 

0.05 

1.01745329 

4.5 

1.43869943 

1.4 

1.26716911 



4.0 

1.42667468 

1.3 

1.25468729 



3.5 

■■ 1.41116081 

1.2 

1.24140742 


. 


6i 


TABLE' 2 


Pj.Q versus A for Pg=2.0 


Pro ^ Pro • ^ ^ro 


30.0 

1.99566531 

25.0 

1.99378314 

20.0 

1.99075733 

15.0 

1.98312290 

13.0 

1.97778900 

10.0 

1.96361072 

9.5 

1.95999964 

9.0 

1.95584057 

8.5 

1.95102074 

8.0 . 

1.94539802 

7.5 

i . 9 3879 165 

7.0 

1.9309696 i 

6.5 

1.92163123 

6.0 

1.91038294 

5.5 

l ! 89670441 

5.0 

1.87990072 

4,5 

1.85903477 

4.0 

1.83283118 

3.5 

1.79954054 


3.0 

1.75674842 

2.0 

1.74677559 

2. -8 

1.73626088 

2.7 

1.72516999 

2.6 

1.71346631 

2.5 

1.70111082 

2.4 

1.68806182 

2.3 

1.67427489 

2.2 

1,65970261 

2.1 

1.64429452 

2.0' , 

1.62799679 

1.9 ' 

1.61075221 

1.8 

1.59249983 

1.7 • 

1.57317490 

1.6' 

1.55270870 

1.5 

1.53102834 

1.4 

1.50805668 

1.3 

1.48371223 

1.2 

1.45790917 


1.1 ■ 

1.43055740 

1.0 

1.40156275 

0.9 

1.37082749 

0.8 

1.33825096 

0.7 

1.30373068 

0.6 

1.26716417 

0.5 

1.22845150 

0.4 

1.18915376 

0.3 

1.14479911 

0.2 

1.09865524 

0.1 

1.05048639 

0.09 

1.04554746 

0.08 

1.04058509 

0.07 

1.03559894 

0.06 

1.03058867 

0.05 

1.02555392 
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TABLE 3 


versus A for ^ “4.0 


A 

J’ro 

A 

^ro 

A 

^ro ■ 

30.0 

3.98647532 

3.0 

3.28185204 

1.1 

2.34942020 

25.0 

'3.98061584 

2.9, 

, 3.25340634 

1.0 

2.26446465 

20.0 

3.96997054 

2.8 

3.22344931, 

0.9 

2.17378415 

15.0 ' 

3.94757229 

2.7' 

3.19188359 

0.8 

2.07670250 

13.0 

3.93112689 

2.6 

3.15860465 

- 0.7 

1.97448896 

10.0 

, 3.88767972 

2.5 

3.12349993 

0.6 

1.86668690 

9.5 

3.87667309 

2.4 ‘ 

3.08644813 

0.5 

1.75226519 

•9.0 

3.86402442 

2.3 

3.04731832 

0.4 

1.62983025 

8.5 

3.84940305 

2.2 

3.00596910 

0.3 

1.49741777 

8.0 

3.83239442 

2.1 

2.96224735 

0.2 

1.35209960 

7.5 

3.81247442 

2.0 

2.91598718 

0.1 

1.18915376 

7.0 

3.78892463 

1.9 

2.86700832 

0.09 

1.17161342 

6.5 

3.76103487 

1.8 

2.81511469 

0.08 

1.15380645 

6.0 

3.72753779 

' 1.7 ■ 

2.76009242 

0.07 

1.13572033 

5.5 

3.68701782 

1.6 

2.70170777 

0.06 

1.11734149 

5.0 

3.63753356 

1.5 

2.63977452 

0.05 

1.09865524 

4.5 

3.57648709 

1.4 

2.57380115 



4.0 

2.50036569 

1.3 

2.50368741 



, 3.5 

3.40436783 

1.2 

2.42902044 
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- TABLE 4 


versus A for Pg=6,0 


A • 

J^ro 

A ■ 

^ro 

A 

^ro 

30.0 

5.97818805 

3.0 

4.85050528 

1.1 

3.35083261 

25.0 

5.96874159 

2.9 

4.80510675 

1.0 

3.21184389. 

20.0 

. 5.95158483 ' 

2.8 

4.75729023 

0.9 

3.06353634 

15.0 

5.91550730 

2.7 

4.70689739 

0.8 

2.90750063 

13.0 . 

5. .88903580 

2. >6 

4.65375768 

0,7 

2.74260186 

10.0 • 

5. 81^6754 

,'2.5 

4.59768678 

0.6 

2.56713263 

9.5 

5.80148207 

2.4 

4.53848518 

0.5 

2.37875492 

• 9.0 ^ 
•,.8.5,. 

5,78116498* 

,, 

, 4.47593656 

0.4 

2,17411590 

5.75768803; 

‘ 2'.2 ' 

'4.40'980585 

0.3 

1.94809778 

8.0’ 

. 5.73038937 

2..1, 

4.33983717 

0,2 

1.69215542 

7.5 

5.69843286 

2.0 

4.26575101 

0.1 

1.38985430 

7.0 

5.66075303 

1.9 

4.18724149 

0.09 

1.35592238 

6.5 

5.61597981 

1.8 

4.10397259 

0.08 

•1-.-32U1922 

6.0 

5,56233483 

1.7 

4.01557400 

0.07 

1.28537407 

5.5 

5.49748709 

1.6 

3.92163593 

0.06 

1.24860602 

5.0 

5.41835022 

1.5 

3.82170298 - 

0.05 

1.21072189 

4.5 

5.32079442 

1.4 

3.71526662 



4.0 

5.19923231 

1.3 

3.60175480 



3.5 

5.04601154 

1.2 

3.48052188 
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TABLE 5 


IJ.Q versus A for =10.0 


A , A Pro- A % 


30.0 

9,. 96230993 

25.0 

9.94598984 

20.0 

9.91635353 

15.0 

9.85405193 

13.0 

9.80835372 

10.0 

9.68779554 

9.5 

9.65729118 

9.0 

9.62225356 

8.5 

9.58177395 

8.0 

9.53471429 

r.5 

9.47963733 

7.0 

9.41471222 

6.5 

9.33758552 

6.0 

9.24520336 

5.5 

9.13356413 

5.0 

8.99736908 

4.5 

8.82952800 

4.0 

8.62044192 

3.5 

8.3'5694816 


3.0 

8.02071861 

2.9 

7.94762822 

2.8 

7.86036927 

2.7 

7.77366594 

2.6 

7.68222058 

2.5 

7.58571070 

2.4 

7.48378688 

2.3 

7.37606757 

2.2 

7.26213835 

2.1 

7.14154576 

2.0 

7.01379247 

1.9 

6.87833134 

1.8 

6.73455792 

1.7 

6.58180171 

1.6 

6.41931492 

1.5 

6.24625905 

1.4 

6.06168748 

1.3 

5.86452363 

1.2 

5.65353259 


1.1 

5.42728355 

1.0 

5.18463272 

0.9 

4.92872952 

0.8 

4.65879094 

0.7 

4.37221815 

0.6 

4.06549504 

0.5 

3.73365882 

0.4 

3.36929763 

0.3 

2.96042648 

0.2 

2.48517268 

0.1 

1.89421267 

0.09 

1.82461982 

0.08 

1.75226520 

0.07 

1.67679133 

0.06 

1.59775623 

0.05 

1.51460253 
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INTRODUCTION 

One of the main difficulties in the use of analytical or computer solu^ 
tions of inherently compensated gas bearings to predict dynamic performance 
is the apparent disagreement with experimental results [1][2], In the 
view of this writer the crux of the problem is that the solutions for 
capacity, stiffness, and damping depend upon the discharge coefficient for 
the inherent orifice as the flow enters the bearing surface from the inlet 
hole. Limited experimental values for the discharge coefficient can be 
found [3], but they are not entirely satisfactory. Furthermore, there is 
no knowledge of the dependence of discharge coefficients on the frequency 
of the disturbance. 

The purpose of this initial study on inherently compensated discharge 
coefficients is to determine experimentally the velocity profiles near 
the feed region of the inlet. Eventually dynamic studies would contribute 
the main effort. This work is only in the beginning stage, and a brief 
outline of the progress to date follows. 


EXPERIMENTAL PROCEDURE 

The characteristic film thickness of gas bearings is at most a few 
thousandths of an inch. Thus, flow visualization studies require a scaling 
based upon the Reynolds number [4], A plexiglass model of the feed inlet 
was constructed with a slot gap adjustable from 1/8 in. to 3/8 in. and a 
width of 2 in.. Figure 1, A water circulating system was built into the 
set up with a pump capacity sufficient to achieve Reynolds numbers on the 


order of 6000. 
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Birefringence Effect 

For many years photoelasticity has been a. very useful method for 
experimental stress analysis in solids. This method is based on the fact 
that certain substances become temporarily doubly refracting when they 
are stressed. An analogous phenomenon in liquids is known as flow birefrin- 
gence. The photoviscous analysis of the birefringence patterns to determine 
velocity distributions in steady, laminar flow is widely known [5]. 

I m > 

When an aqueous solution of commercial organic dye. Milling Yellow NG5, 
is caused to flow through a transparent viewing channel and subjected to 
polarized light, the fluid becomes double refracting and produces visible 
interference patterns which are seen as alternating dark and light bands . 

The dark bands, iso chromatics, represent the loci of points where the 
magnitude of the shearing stresses are the same. The velocity profiles 
are then obtained from integration of the shear stress profiles. 

The optical system used in the experiments was readily available in 
the experimental stress lab at Mississippi State University. The basic 
system consists of a sodium-vapor light source, two collimating lenses, 
two polarizing plates, two quarter wave plates, and an optical bench. 


Particulate Effect 

An alternate source of flow visualization consisted of particles j 
suspended in water. Motion of these particles was recorded in the following 
manner. A collimated light beam with a width of approximately 1/4 of an 
inch was projected vertically through the flow from above. Figure 1. When 
viewed from a horizontal position, the motion, of these particles was visible 
.in this plane of light. Photographic records of particle motion were made. 

, Two particulates were used: (1) small hollow glass spheres approximately 

i • t * 

’(. 001 inch in size, (’2) a certain 'dishwater detergent which apparently has 

^ it ± T * K i * 




FIG. 1. Visualization of Flow Over a Corner. 


o 
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suspended extremely small abrasives. Trial and error procedures lead to 
the use of 1/100 gram per gallon of the former and 2 drops per gallon of 
the latter. 

Hot Wire Anemometer 

In place of the aqueous solutions above a direct attack on the velocity 
profiles can be obtained from hot wire measurements using air. This proce- 
dure is so standard that no futher summary is given. The biggest drawback 
to the use of a hot wire is the subsequent interference of the flow field 
by the probe itself — an unknown effect. 


RESULTS 


Birefringence Effect 

The original plan was to, make use of a bir-efringent fluid since the 

r ' t 

flow field is unobstructed and the field can be obtained throughout with 
one test run. Unfortunately the literature on birefringence was not explicit 
in explaining several difficulties when working with Milling Yellow. This 
approach was abandoned for the following reasons: 

(1) Very little is known about the rheological properties of Milling 
Yellow, particularly the relationship of birefringence with 

f 

the shear stress which is non-Newtonian; 

(2) concentrations as narrow as 1 1/4 - 1 1/2 per cent are necessary 
and these concentrations show a highly nonlinear relationship 
with the amount of birefringence; 

(3) birefringence is very sensitive to temperature with as little 
as one degree variation introducing a sizable error. Pumping 
losses can account for several degrees change in the solution 


temperature; 
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(4) the planned dynamic studies would undoubtably result in a 
total effort to obtain the rheological properties of the 
Milling Yellow and very little realizable data on the fluid 
mechanics of the problem. 

Particulate Effect 

After a concerted effort with the particulate solution approach to 
flow visualization this phase was abandoned for the following reasons: 

(1) Although the flow patterns were observable to the naked eye, 

it proved to be very difficult to photograph the patterns, 

and professional photographers at the University offered 

* 

little hope; 

(2) the experience also suggested that the photographs, if 
obtainable, could be used for a qualitative picture of the 
flow field but of marginal value determining the velocity 
profiles . 

Hot Wire Anemometer 

The- hot wire anemometer appears to be the only logical choice for 

establishing the flow field accurately. A test run was conducted by simply 

inserting a miniature probe into the slot opening and tracing the flow 

field from the leading edge of the step. Pressure taps were placed at 

several stations from the leading edge so that the discharge coefficient 

could be found directly. Tenative results were compared with the analytical 

work of Hagerup [4]. Some of the findings are as follows. 

A well defined separation bubble exists and its length is approximately 

two film thicknesses from the step edge. Figure 2.' This value corresponds 
; * 

to a theoretical Reynolds number of 500, Figure 3. The calculated value 
for the test run was approximately 3600 based upon an average exit velocity. 
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Calculated separation atreaialines for 
several Reynolds numbers 



Inherent restrictor pressure loss coefficient 
as function of Reynolds number 


FIG. 3. Theoretical Results from Hagerup [41 
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Hagerup [4] defines the loss coefficient as 
Ap = "I kpU^ 


which has a corresponding discharge coefficient 


c = [— ?•■ -1^^^ 
d ^1+k^ 


The theoretical discharge coefficient for R = 500 is c, = 0.88. The calculated 

a 

values from the experiment were ^ 0.88 for a Re 3 rnolds niimber range 2400- 
4000. 

The results would agree with the theoretical calculations provided the 
measured Reynolds number was smaller by a factor of 6-8. All velocity ^ 

measurements were referenced to a ”Magnehelic” pressure gauge with a direct 
read out in inches of water. There is insufficient data to draw meaningful 
conclusions from the experiments . It is possible that the probe itself is 
interf erring with the flow pattern. Future experiments will be conducted 
with a larger model of the inlet corner. 



76 


REFERENCES 


1. Pink, E, G., "Investigations into Design Methods for Externally Pres- 
surized Gas Journal Bearings", 6th" Inter. Gas Bearing S3nn. , March 27-29, 
1974, Univ. of Southampton. 

2. Cunningham, R. E., et.al., "Design of a Squeeze Film Damper for a Multi- 
Mass Flexible Rotor", J. of Eng. for Ind., ASME Trans., Nov. 1975, 

p. 1383-1390. 

3. Vohr, J. H., "An Experimental Study of Flow Phenomena in the Feeding 
Region of an Externally Pressurized Gas Bearing", MTI-65TR47, Mechanical 
Technology, Inc., 1965. 

4. Hagerup, H. J., "On the Fluid Mechanics of the Inherent Restrictor", 

ASME paper 74-FE-26, 1974. 

5. Wayland, H., "Streaming Birefringence as a Rheological Research Tool", 

J. of Polymer Science, Part C, No. 5, 1964, p. 11-36. 




